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DIGITAL  CONTROL  OF  FLIGHT  SIMULATOR  MOTION  BASE 


by 

J.  SANDOR 


\ 

Digital  control  strategies  for  a  nonlinear  motion  base 
actuator  are  considered  and  a  compound  linear/nonlinear  algorithm  is 
derived  for  velocity  tracking  under  a  range  of  load  conditions  for 
representative  motion  demands.  Discrete  frequency  domain  methods  are 
employed  to  synthesize  the  linear  components  of  the  control  law  and 
to  establish  the  stability  of  the  closed  loop  system  via  describing 
function  analysis.  The  nonlinear  compensation,  providing  attenuation 
of  Coulomb  friction  effects,  is  realised  through  a  velocity  dependent 
gain  term,  which  does  not  significantly  affect  the  satisfactory  large 
amplitude  system  response.  Overall  closed  loop  system  performance  is 
validated  through  computer  simulation. 
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1.  INTRODUCTION 


Motion  cue  generation  is  a  significant  factor  is  any  manned 
flight  simulation  facility.  Such  a  facility  is  being  developed  at  the 
Aeronautical  Research  Laboratories  in  support  of  research  in  synthetic 
training  technology  and  aircraft  system  studies.  A  key  element  of  the 
development  work  is  the  formulation  of  (digital)  control  strategies 
for  motion  base  systems  to  synthesize  position,  velocity  and  acceleration 
in  a  manner  that  replicates,  most  effectively,  the  cueing  environment, 
subject  to  the  overall  excursion  and  rate  constraints  intrinsic  to  the 
motion  base  actuators.  In  support  of  these  aims  investigations  are 
being  made  into  the  synthesis  of  tunable  control  laws  for  velocity 
control  of  actuators  under  a  range  of  load  and  reference  input 
conditions  in  the  presence  of  intrinsic  nonlinearities  and  time  delays. 

This  report  describes  the  derivation  of  a  control  algorithm 
for  the  electric  servomotor  drives  of  a  4  degrees-of-freedom  synergistic 
simulator  motion  platform.  The  performance  criteria  are  specified  in 
terms  of  the  error  between  the  reference  velocity  input  and  the  motor 
velocity  response  profile  for  a  range  of  inertial  and  spring  restoring 
loads  and  subject  to  maximum  allowable  response  time  delay.  (Here 
reference  input  velocity  refers  to  the  motion  base  drive  signal  and 
includes  washout) .  In  addition  to  the  preceding  requirements ,  control 
law  complexity  is  constrained  to  be  such  that  a  four  channel  controller 
can  be  realized  on  a  single  DEC  LSI-11/23  microcomputer  operating  at 
50  Hz  sampling  rate. 

The  control  law  presented  comprises  a  linear  proportional  plus 
integral  (PI)  term  with  compensation  to  meet  performance  and  tunability 
requirements  and  a  nonlinear  term  to  attenuate  the  transient  effects  of 
Coulomb  friction.  This  control  regime  is  derived  in  three  steps: 
firstly  a  discrete  frequency  domain  approach,  based  on  a  z-transform 
model  of  the  linearized  system, is  adopted  for  the  derivation  of  the 
linear  terms  of  the  controller  such  that  satisfactory  stability  margins 
are  achieved  for  the  linearized  system  under  the  specified  load  conditions 
and  loop  time  delays.  Secondly,  the  stability  of  the  nonlinear  closed 
loop  system  is  established  by  using  harmonic  linearization  for  mathe¬ 
matically  characterising  Coulomb  friction  and  the  application  of  describing 
function  analysis.  Finally,  based  on  the  simulation  of  the  linearly 
controlled  system,  a  smooth  nonlinear  control  term  is  synthesised  for 
the  attenuation  of  friction  effects  which  are  shown  to  be  particularly 
significant  for  low  amplitude  reference  inputs.  The  nonlinear  control 
term  effects  an  increase  in  the  system  loop  gain  in  regimes  where  Coulomb 
friction  becomes  dominant,  without  adversely  affecting  the  satisfactory 
large  aaplitude  response  of  the  linearly  controlled  system. 

In  the  present  case,  given  the  controller  ooeqplexity  constraints 
and  the  requirement  for  tunability,  the  frequency  domain  synthesis  method 
is  preferred  to  that  of  either  optimal  control  theory  or  deadbeat  controller 
design.  The  former  is  based  on  the  minimization  of  a  quadratic  penalty 
function,  which  here  would  require  the  solution  of  a  computationally  costly 
two  point  boundary  value  problem  [9],  whilst  the  realization  of  the 
resulting  control  algorithm  requires  estimation  of  all  the  system  states. 


Deadbeat  control  algorithms,  on  the  other  hand  require  a  priori  knowledge 
of  system  inputs,  are  not  readily  tunable,  and  may  require  excessive 
control  effort  to  achieve  deadbeat  response  [4].  Zn  contrast,  frequency 
domain  methods  offer  readily  computable,  tunable,  time-invariant  control 
laws,  implementable  in  recursive  form  and  giving  satisfactory  response 
for  a  range  of  system  inputs  and  loads.  Their  principal  limitations, 
apart  from  non-optimality,  is  that  for  nonlinear  systems  closed  loop 
stability  is  not  a  priori  guaranteed.  However,  for  the  system  under 
consideration  here,  where  the  linear  parts  of  the  system  have  low  pass 
frequency  characteristics,  this  difficulty  can  be  overcome  by  application 
of  describing  function  analysis. 


2.  SYSTEM  MODEL  AND  PERFORMANCE  SPECIFICATION 
2.1  Motion  Base  Description 

The  four  degrees-of-freedom  simulator  motion  base  under 
consideration  comprises  a  metal,  two-sum,  side  by  side,  cockpit  mounted 
on  a  pneumatically  balanced  suspension  system  with  motion  actuation 
derived  from  4  direct  current  printed  circuit  motors  via  screw  jacks 
attached  to  the  oockpit.  Two  motor /screw jack  actuators  are  attached  to 
the  rear  of  the  cockpit  whilst  another  pair,  inclined  to  the  vertical, 
jure  attached  through  a  "knuckle"  arrangement  to  the  front  of  the  cockpit 
as  depicted  in  Figure  2.1.  The  operation  of  the  motion  base  is  synergistic 
that  is,  it  requires  a  combination  of  actuator  extensions  to  achieve 
motion  in  any  one  degree  of  freedom.  This  "coupling"  effect  results  in 
v jurying  inertial  and  spring  restoring  loads  on  each  motor.  For  example, 
heave  motion  is  realized  with  all  four  actuators  driving  the  load,  which 
includes  the  spring  restoring  force  of  the  suspension  system.  On  the 
other  hand,  for  yaw,  the  front  actuators  drive  the  load,  with  the  rear 
actuators  only  tracking,  and  the  suspension  system  providing  practically 
no  restoring  force. 

The  screw  jacks  comprise  a  25  mm  lead  screw  and  a  nut  with  a 
maximum  linear  travel  of  415  nn  for  the  front  jacks  and  590  mo  for  the 
rear  jacks.  The  maximum  oockpit  excursion  about  the  central  position 
for  each  degree  of  freedom  is  - 

heave  :  i  300  mn 
pitch  s  ±25° 
roll  s  ±17° 
yaw  :  ±  12" 

Each  motor  is  a  Yaskawa  printed  circuit  d.c.  motor,  model  UGPMFE  -  16AAB 
rated  at  350W,  47V.  The  weight  of  the  cabin  fully  loaded  is  460  kg. 
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2.2  Motor/Load  Model 


If  a  voltage  V  is  applied  to  a  motor  with  armature  resistance 
r  the  gross  electrical  torque,  T,  is  given  by  - 

A 

T  *  k  I  ;  V  *  E  +  r  I 
T  a 

where  I  is  the  current  in  amps ,  E  is  the  e.m.f.  in  volts  and  k  is  the 
motor  torque  constant.  Noting  that  at  speed  w  the  e.m.f.  constant  k 
is  equal  to  E/w  and  since  electrical  power  available  for  conversion  E 
must  equal  gross  mechanical  output  power  (El  =  Tw)  the  motor  equation 
reduces  to  - 


r 

V  *  k  w  +  T  . 

T  k 

T 

Now,  gross  torque,  T,  is  given  by  - 


.. (2.1) 


T 


V 


+  T,  +  TV 
P  0 


.. (2.2) 


where  k^  is  the  damping  torque  constant,  T  is  the  Coulomb  friction 
torque  and  Tq  is  the  available  output  torque.  Combining  (2.1)  and  (2.2) 
the  speed-torque  characteristic  of  the  motor  is  - 


where 

1/f 


is  motor  regulation  (reduction  in  speed  per  unit  change  in  load  torque) . 
Typically  for  motors  over  300W,k  is  2%  of  (k  a/r  )  (1]  and  hence  f  can 
be  taken  as  -  D  T  a 


k  *  +  k 
T  D 


f  i  k  a/r  . 
T  a 


Now  for  an  inertial  load  J,  and  spring  restoring  load  torque 
k^fl,  where  8  is  the  motor  angle,  the  dynamic  equation  of  the  motor  is  - 


J 


dw 

dt 


fw 


k  6 
s 


Taking  Laplace  transforms,  the  linear  system  transfer  function  is  - 


running  friction,  T  *  0.2  Nm. 


2.3  Control  System  Configuration 

Figure  2.3  shows  the  computer  control  configuration  for  each 
servomotor  at  a  sampling  period  of  At  seconds.  The  control  computer 
sauries  the  reference  input  and  system  output  both  of  which  are  available 
in  digital  form.  The  reference  velocity  input  is  computed  by  the  flight 
simulation  oomputer  whilst  the  control  computer  calculates  the  motor 
velocity  from  the  output  of  the  shaft  position  encoder  which  has  a 
resolution  of  1/600  revolution  of  the  motor  shaft.  It  will  be  assumed 
that  the  encoder  is  linear  with  bandwidth  sufficiently  great  compared 
to  the  system  bandwidth  such  that  its  dynamics  need  not  be  considered 
in  the  subsequent  analysis. 

In  operation  the  control  computer  calculates  the  system  velocity 
error  and  computes  the  control  which  is  to  be  applied  to  the  system,  via 
a  zero  order  hold  circuit,  one  sampling  period  later.  This  computation/ 
control  delay  is  represented  in  the  model  by  the  transfer  junction 
“SAt 

e  .  Although  the  d.c.  servo  amplifier  is  represented  in  the  config¬ 
uration  of  Figure  2.3,  it  will  be  assumed  to  have  a  bandwidth  such  that 
its  analytic  representation  is  just  a  pure  gain  in  the  control  law. 

2.4  Performance  Specification 

The  aim  of  the  control  system  design  is  to  derive  a  control 
algorithm  that  will  maintain  closed  loop  system  stability  and  performance 
under  a  range  of  system  parameters  and  load  conditions.  In  general  the 
specification  of  the  performance  of  a  motion  base  platform  involves  the 
characteristics  of  the  actuators,  platform  dynamics,  the  aircraft  to  be 
simulated  and  the  physiology  of  human  motion  sensory  system.  In  the 
present  case  consideration  is  limited  to  the  individual  actuator  loops 
and  hence  control  system  performance  can  be  adequately  specified  in 
terms  of  accepted  servomechanism  criteria. 

Specifically,  for  a  sampling  period  - 


At  *  0.02  sec. 


and  for  nominal  system  parameter  values  of  - 

J  *  4  x  10  ^  kg  ma 

k  ■  2.5  x  10  ^  Nm/rad 

8 

the  compensated  open  loop  system  is  to  have 

(i)  8db  gain  margin; 

(ii)  45s  (minimum)  phase  margin;  and 

(iii)  closed  loop  system  bandwidth  (3db)  of  5  Hz. 


DIGITAL  CONTROL  CONFIGURATION 
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• 

The  choice  of  sampling  interval  is  primarily  dictated  by  the  computational 
capability  of  the  LSI-11/02  subject  to  the  requirement  that  a  multiplexed 
four  channel  control  system  is  to  be  realized.  The  gain  and  phase  margin 
specifications  reflect  closed  loop  system  stability  with  respect  to 
system  parameter  variations  whilst  closed  loop  bandwidth  is  derived  from 
the  reference  input  bandwidth  which  is  estimated  to  be  no  greater  than  .  • 

2  Hz.  The  bandwidth  requirement  will  ensure  adequate  speed  of  response. 

In  addition  to  the  above  requirements  a  further  measure  of 
the  tracking  accuracy  of  the  closed  loop  system  can  be  specified  as 
follows:  let  the  time  translated  sum  of  squares  velocity  error  be 

defined  as  -  • 


S  (j)  A  I  (  w(kAt)  -  w  (k+j  At)  ) 

W  £  T 

k=0 

and  define  the  "average"  time  delay  of  response  as  -  • 


average  time  delay  =  Lht 


where 


l  =  min  S  (j)  . 

j 


Then  a  "global"  measure  of  distortion  is  given  by  the  rms  quantity  - 


where  NAt  is  the  maximum  of  the  duration  of  the  motor  response  and  the 
reference  input.  Defining  S. ,  the  time  translated  sum  of  squares 
acceleration  error,  analogous  to  S  ,  a  measure  of  acceleration  tracking 
error  is  given  by  - 


R. 


w 


! 

(kAt)3 


and  reflects  the  fidelity  with  which  accelerations,  and  hence  forces, 
can  be  imposed  upon  the  crew  in  the  simulator. 

Since  the  present  study  is  limited  to  actuator  control  rather 
than  total  platform  oontrol ,  bounds  on  rms  tracking  errors  are  given  only 
for  sample  reference  inputs  derived  below. 


_» 


-  9  - 


Specifically  these  bounds  are  - 


<  0.1 


Rw  U)  <  0.5 


Iht  K  60  msec. 


The  latter  bound  is  set  so  as  to  allow  adequate  time  for  reference  input 
computation  and  still  keep  the  total  cue  generation  time  below  100  msec; 
the  delay  time  often  quoted  for  the  onset  of  pilot  induced  oscillation. 
On  the  other  hand  the  bounds  on  rms  tracking  errors  are  only  estimates 
constrained  by  encoder  accuracy;  a  more  expedient  and  elegant  approach, 
which  is  to  be  employed  in  the  design  of  the  overall  control  system  for 
the  motion  platform,  is  to  include  the  tracking  error  in  the  performance 
index  of  an  optimal  control  formulation  of  the  control  law  design. 

Finally,  the  threshold  velocity  error  limit  at  which  the 
actuator  is  to  pull  free  from  static  friction  is  taken  as  0.67  rad/sec. 

2.5  Sample  Reference  Velocity  Inputs 


The  time  profile  of  motion  platform  cues  consist  of  two  phases: 
the  first  is  the  "alerting"  or  "onset"  phase  during  which  the  platform 
motion  produces  pilot  sensed  motion  cues,  whilst  during  the  second  phase 
or  "washout"  the  platform  is  driven  below  subliminal  rates  back  to  the 
neutral  position.  The  choise  of  onset  time,  magnitude  of  a  motion  cue 
and  washout  is  a  compromise  between  maximum  velocity,  acceleration  and 
extension  of  platform  actuators  and  the  subliminal  levels  of  motion 
perception.  That  is,  for  washout  at  below  subliminal  acceleration  and 
velocity  rates,  a  large  onset  cue  requires  long  actuator  travel. 

For  the  present  purposes  of  actuator  controller  design  two 
reference  velocity  inputs  have  been  derived  based  on  an  onset  time  of 
0.25  sec,  maximum  actuator  travel  equal  to  75  percent  of  maximum  extension, 
with  the  subliminal  acceleration  and  velocity  levels  quoted  in  [2]  taken 
as  being  indicative  of  thresholds  of  motion  sensation.  The  reference 
inputs,  w  and  w  ,  are  shown  in  Figure  2.4  with  w  ,  corresponding  to 


inputs,  w  and  w  ,  are  shown  in  Figure  2.4  with  w  corresponding  to 
motion  inrlthe  vertical  axis  and  w  to  motion  in  one  of  the  rotational 
degrees  of  freedom.  £ 


of  the  rotational 


CONTROL  COMPENSATOR  SYNTHESIS 


3.1  Basic  Considerations 


As  a  starting  point  for  digital  controller  design  a  proportional 
plus  integral  (PI)  control  law  of  the  form  - 


V(s)  -  k  (  1  +  —  )  <n_(s)  -  fl(s)) 
P  s  R 


.. (3.1) 


is  chosen  for  the  continuous  time  representation  of  the  motor/load 
system  excluding  the  effects  of  sampling  and  friction  as  depicted  in 
Figure  3.1. 


FIG.  3.1  CONTINUOUS  -  TIME  MOTOR/LOAD  MODEL 


Derivative  control,  which  would  add  an  extra  degree  of  freedom,  is  not 
considered  since  neither  a  sufficiently  accurate  estimate  of  acceleration 
is  obtainable  from  shaft  encoder  measurements  nor  is  it  possible,  in 
the  present  scheme,  to  measure  armature  current.  The  proportional  and 

integral  constants  k  and  k.  are  based  on  - 

P  i 

(i)  pull  out  torque  of  the  motor  to  overcome 
stiction;  and 

(ii)  dominant  closed  loop  pole  placement  that 

is  relatively  insensitive  to  spring  rate, 

k  ,  including  the  case  k  *  0. 
s  s 

For  nominal  values  of  system  parameters  - 

J  *  4  x  10  2  kg  m3 

f  *  1.82  x  10  2  Nm/rad/sec 

k  *  2.5  x  10  2  Nm/rad 

s 

the  proportional  gain,  k_,  is  selected  such  that  the  motor  torque  is 
sufficient  to  overcome  the  measured  stiction  torque  at  the  velocity  error 
threshold  specified  in  Section  2.4,  i.e. 

S 

*P  ~  <vr  '  wl  threshold  *  *tiction  tot^e 
a 


13  - 


and  hence 


k  •=  2.5  volt/volt. 

P 


One  the  other  hand  root  locus  analysis  for  the  two  systems  k  >0  and 


k  -  2.5 
tiiat  - 


x  10-2  (Figure  3.2) 


suggests  that  k^ 


should  be  chosen  such 


k.  >  6 

l 


if  the  closed  loop  characteristics  of  the  system  are  to  be  similar  for 
the  two  spring  rate  conditions.  For  the  present  design  k.  is  chosen 
as  - 


k.  *  15. 


.  • 


.  • 


lb  gauge  the  effects  on  closed  loop  stability  arising  from 
digital  realization  of  the  PI  control  law,  computational  time  delay  and 
sampling  have  to  be  taken  into  consideration.  For  this  purpose  consider 
the  model  of  the  closed  loop  system  depicted  in  Figure  3.3  where 
exp(-sT)  and  (1  -  exp(-sT))/sT  are  three  transfer  functions  for  computer 
time  delay  and  the  sampler  plus  zero-order  hold  circuit  [6].  Defining 
the  open  loop  transfer  function  as  - 

k  kT  (1  -  exp(-sT) )  Z' 

To(s)  “  ~T~  exp(_sT)  - ^ -  (^1  + 

then  closed  loop  stability  can  be  assessed  from  the  Nyquist  frequency 
of  T^ts)  13].  For  the  values  of  system  and  control  parameters  the 
Nyqurst  plot  of  T  (ju>)  is  shown  in  Figure  3.4.  Recalling  that  stability 
of  the  closed  loop  system  is  assured  if  the  graph  of  T  (ju>)  does  not 
make  a  counter-clockwise  encirclement  of  the  -1+jo  point,  it  cam  be 
seen  from  Figure  3.4  that  this  condition  is  not  satisfied,  amd  hence 
the  closed  loop  system  of  Figure  3.3  will  be  unstable  amd  additional 
compensation  is  required.  The  instability  is  due  to  the  significant 
phase  shift  introduced  by  the  time  delay  and  saunpler/zero-order  hold 
circuit  as  shown  in  Figure  3.5. 

3.2  Digital  Control  Law 

Since  the  controller  is  to  be  realized  digitally,  the  most 
convenient  system  description  is  in  terms  of  z-transforms  [3], [4]. 
Adopting  the  convention  that  if  g(s)  is  the  Laplace  transform  of  g(t) 
then  G(z)  is  its  z-transform  (with  z  ■  esT)  the  closed  loop  digitally 
controlled  system  is  depicted  in  Figure  3.6  below.  Here  account  is 
taken  of  the  fact  that  motor  speed  is  derived  from  sampled  position 
measurements  by  simple  differencing  and  reference  input  position  is 
computed  by  rectangular  integration  of  reference  velocity  input. 


Hi') _ ! _ 

s  J  Js2  +  fs  +  k 


FIG.  3.3  SYSTEM  MODEL  INCLUDING  TIME  DELAY,  SAMPLING 
AND  ZERO  ORDER  HOLD 
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Let  G(s)  be  the  composite  transfer  function  of  the  zero-order 
hold,  motor/load  and  integrator  with  output  being  taken  as  (the  physically 
instrumented)  motor  position.  Then  - 


G(s) 


rJ  sa  +  2  5  w  p  +  u)a*s" 
a  n  n 


_1  1  -  exp(-sT) 

s 


where 


/ 


k  /J 
s 


C  = 


2  A  j1 

s 


.. (3.3) 


Then  a  straight  forward  but  messy  calculation  [3]  gives  the  z-transform 
of  G(s)  as  - 


(i)  for  k  y  0,  5  <  1 
s 


G(z) 


T  f  ,  -6T 

- (  i  -  iil® _ 

r  u  aJ  V  ,  o 

an  V  za  -  2: 


(cos(yT)  -  (6/y) 


-6T 

2ze  cos (yT) 


y)sin(yT)))  (z-P^N 

♦  „-2{T  J 


. (3.4a) 


and 


6  =  £w 


n 


y  =  w  /I  -  5 2 '  ; 
n 


(ii)  for  k  f  0,  c  >  1 
s 


G(z) 


3—  A  »  HjlIL  (  _i _ « _ ^ 

WV  Vz-e-5T 


. . (3.4b) 


and 


<5  =  ujr  (?  +  ) 


y  -  uin  (;  -  ) 


SAMPLED  -  DATA  MODEL  OF  SYSTEM 


E 


for  k 


1  +  Y)  +  (1 


r  Jo3  z3  -  (1  +  y)z  +  y 
a 


Computing  motor  velocity,  w,  by  simple  differencing  of  motor  position 
and  confuting  reference  input  position,  9  ,  by  rectangular  integration, 
the  z-transform  of  the  control  law,  incluSing  computation  delay  is  - 


L(.)  =  ^  +  ki  > 

z  zT 


.. (3.5) 


Taking  the  form  of  compensation  to  be  a  second  order  section  with 
transfer  function  - 


Cb„z3  +  b_z  +  b,  N 

•  •:) 


.  (3.6) 


then  the  closed  loop  system  transfer  function  is  given  by  - 


fcl(z) 


F  (z) 
o 

1  +  F  (z) 
o 


. (3.7) 


F  (z)  *  L (z)  D(z)  G(z). 

o 


. (3.8) 


The  choice  of  compensator  parameters  is  dictated  by  examination 
of  closed  loop  system  stability.  This  is  accomplished  by  evaluating  the 
open  loop  transfer  function  F  (z)  along  the  unit  circle  |z|  =1. 
Equivalently,  via  the  bilinear  transformation  - 


r  +  1 

F^T  •  r  ■  j“r 


the  stability  of  the  closed  loop  system  is  determined  by  the  Nyquist 

plot  of  F  (jw  ), 
o  J  x 


/'ju  +  IN 
£  Fo  (  ) 
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• 

as  >*)  varies  from  -®  to  0  [4]  ,  where  u>  can  be  viewed  as  a  "generalized 
frequency".  Analogously  to  the  continuous  time  case,  counter-clockwise 
encirclement  of  the  -1  +  jo  point  by  F  (j-'  )  implies  instability  of  the 
closed  loop  system. 

An  interactive  computer  program  with  graphics  capability  was  .  • 

utilized  for  stability  analysis  and  compensator  design.  Firstly,  to 
counter  the  phase  lag  introduced  by  sampling  and  computation  delay, 
lead  compensation  was  considered  with  parameters  - 

a  =  b  =  0,  a  =  1.0  b  =  0.4  a  =  b  *  -0.33 

XX  4  |  4  f  J 

-3  -2 

for  J  =  4x10  ,  k  =  2.5x10  .  The  Nyquist  plot  of  this  design  is  shown 

in  Figure  3.7(a).S  Clearly  the  closed  loop  system  is  stable  but  the  gain 
margin  is  only  1.6  db  and  there  will  be  a  resonance  at  u  =2.  Moreover 
for  small  inertia  (J  =  2xl0-2)  lead  compensation  alone  will  not  stabilize 
the  system  as  shown  in  Figure  3.7(b).  Hence  increasing  the  order  of  the 
compensator  is  required  to  give  additional  degrees  of  freedom  to  maintain  .  • 

phase  lead  and  reduce  gain  at  the  phase  crossover  frequency.  A  design 
choice  for  D(z)  of 

a  =  12  a  =  -10 

1  2 

b1  -  2.2  b2  -  -2.5 

gives  the  required  improvement  in  gain  margin  and  stability  to  low  inertial 

load  as  shown  in  Figure  3.8.  Further  the  stability  of  the  closed  loop 

system  in  the  presence  of  variations  in  load  inertia  and  spring  constant 

is  established  by  the  Nyquist  plots  of  Figure  3.9.  The  closed  loop 

frequency  response  for  the  control/compensator  design  and  J  *  4x10”^, 

k  «  2.5xl0-2  is  obtained  by  plotting  magnitude  and  phase  of 
s 


fcl  <z) 


j2ufT 


as  in  Figure  3.10,  with  3db  bandwidth  being  7  Hz  and  45°  phase  shift  at 
2  Hz. 


3.3  System  Stability  in  the  Presence  of  Static  and  Sliding  Friction 

To  complete  the  stability  analysis  of  the  closed  loop  system 
consideration  has  to  be  given  to  the  presence  of  nonlinearities.  In  the 
present  case  the  dominant  nonlinearity  in  the  loop  is  static  and  sliding 
friction  (since  the  screw  jacks  are  typically  under  load,  backlash  can 
be  neglected) .  The  method  of  analysis  to  be  employed  to  determine  closed 
loop  stability  of  the  resultant  nonlinear  system  is  based  on  describing 
functions  (OF)  or  harmonic  linearization. 

The  salient  features  of  this  approach  can  be  described  by 
consideration  of  the  nonlinear  continuous-time  system  of  Figure  3.11. 

The  method  is  based  on  the  assumption  that  if  a  limit  cycle  occurs  in 


FIG.  3.10(a)  CLOSED  LOOP  FREQUENCY  RESPONSE  -  GAIN 


FIG.  3.10(b)  CLOSED  LOOP  FREQUENCY  RESPONSE  -  PHASE 


the  closed  loop  system,  then  as  most  control  systems  have  low  pass 
frequency  characteristics,  the  input  to  the  nonlinearity  will  be 
approximately  sinusoid  al[5].  The  OF  N(x)  is  the  ratio  between  the 
input  (x)  and  output  amplitude  (xN(x))  of  the  sinusoidally  excited 
nonlinearity.  If  the  nonlinearity  is  double  valued  then  the  OF  is 
complex.  Stability  of  the  closed  loop  system  is  then  evaluated  by 
plotting  G  (jw)G2(jw)  and  1/N(x)  on  a  Nyquist  diagram  and  applying 
the  Nyquist  criterion  with  respect  to  the  locus  -1/N(x)  instead  of 
the  point  -1+jo.  That  is, 16}  if  - 

(i)  the  G  ( jii))G_  ( jijj  curve  encloses  the  -1/N(x)  locus 
the  closed  loop  system  is  unstable; 

(ii)  the  G  (joj)G  (jo))  curve  does  not  enclose  the  -1/N(x) 
locus'1  then  the  closed  loop  system  is  stable;  and 

(iii)  the  G  (jui)G  (juj)  curve  intersects  -1/N(x)  locus  then 

a  limit  cycle,  with  amplitude  and  frequency  determined 
by  the  point  of  intersection,  may  arise.  The  stability 
of  the  limit  cycle  is  ascertained  by  variational 
techniques  15] . 

In  the  case  of  digital  nonlinear  systems  describing  function 
analysis  can  be  carried  out  in  one  of  two  ways:  firstly  by  application 
of  the  discrete  describing  function  method  developed  in  [4,  Chp.  12] , 
or  secondly  by  approximating  the  digital  elements  of  the  control  loop 
by  frequency  response  equivalents  and  applying  conventional  describing 
function  stability  criteria.  The  first  approach  is  specific  to  systems 
with  input  nonlinearities  (such  as  relays)  and  requires  the  examination 
of  possible  limit  cycles  at  each  integral  subharmonic  of  the  sampling 
frequency.  On  the  other  hand  the  alternative  approach  is  only  valid 
if  the  digital  elements  of  the  loop  are  separated  from  the  nonlinearity 
by  continuous-time  elements  having  low  pass  characteristics.  In  the 
present  case,  noting  that  the  nonlinearity  in  question  is  not  at  the 
motor  input  and  that  a  zero  order  hold  is  a  low  pass  element,  the  second 
of  the  two  approaches  is  adopted. 

Consider  the  approximate  model  of  the  closed  loop  system 
depicted  in  Figure  3.12  where  - 


D(jw)  -  L(z)  D(z) 

joiT 

z  ■  e 

is  the  frequency  response  of  the  digital  controller  including  the 
computational  delay. 


FIG.  3.12  FREQUENCY  RESPONSE  CHARACTERISATION  OF 
CLOSED  LOOP  SYSTEM 


Assuming  the  Coulomb  friction  characteristic  of  Figure  3.13,  the 
describing  function  for  friction  acting  on  a  purely  inertial  load  has 
been  derived  in  [8]  and  is  plotted  as  1/N(A)  in  Figure  3.14.  Here  the 
describing  function  parameter  is  - 


X 


VT 


a 


with  the  applied  torque  being  T  sinut  and  the  loci  terminate  at 
A  =  Tj/Ts  since  X  >  T  /T  implies  T  <  T  .  Rearranging  the  block 
diagram  of  Figure  3.1?  such  that  friction  acts  on  an  inertial  load 
leads  to  the  system  representation  of  Figure  3.15(a)  and  further 
rearrangement  results  in  the  system  of  Figure  3.15(b),  where 


G(jw) 


D(jw) 


(1  -  e~iaT) 
ju)T 


+ 


jf*)  • 


Application  of  the  Nyquist  stability  criterion  to  the  plot 
of  1/N(A)  and  b(j^)  *  -k,Vr  G(ju)/Juia  on  the  complex  plane  (Figure 
3.16)  readily  verifies  tfiatathe  closed  loop  system  will  be  stable 
for  the  expected  range  of  system  and  Coulomb  friction  parameters. 


Jl 


FIG.  3.15  SYSTEM  TRANSFORMATION  FOR  DESCRIBING 
FUNCTION  ANALYSIS 
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FIG.  3.16(a)  DESCRIBING  FUNCTION  AND  FREQUENCY  LOCI 
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FIG.  3.16(b) 
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4_. _ SIMULATION  OF  CLOSED  LOOP  SYSTEM 

The  efficacy  of  the  digital  controller  design  in  meeting  the 
required  closed  loop  system  performance  criteria  is  confirmed  by  digital 
simulation.  The  motor  model.  Coulomb  friction  and  position  encoder 
characteristics  are  represented  as  a  set  of  first  order  differential 
equations  with  algebraic  constraints,  and  integrated  using  a  4th  order 
Runge-Kutta  integration  routine  with  a  step  size  of  O.Ol&t  on  a 
DECsystem-10  computer.  Included  in  this  representation  is  the  time 
delay  between  measurement  of  system  states  and  reference  input,  and 
application  of  the  control  input.  The  digital  compensator  is  simulated 
in  recursive  form,  operating  at  a  simulated  sampling  period  of  it  seconds 
with  motor  position  measurements  quantized  to  600  counts  per  revolution 
to  emulate  the  position  encoder.  In  this  way  the  motor  appears  as  a 
continuous  system  to  the  digital  controller  with  system  nonlinearities 
being  faithfully  represented. 

Simulations  were  conducted  over  a  range  system  parameters 
for  the  class  of  reference  inputs  specified  in  Section  2.5.  The  results 
of  these  are  depicted  in  Figure  A1  and  the  computed  performance  measures 
are  tabulated  in  Table  A1  of  the  Appendix.  It  can  be  readily  seen  that 
for  a  wide  range  of  system  parameters  the  closed  loop  system  exhibits 
good  velocity  tracking  and  the  performance  measures  are  within  the 
required  bounds.  However,  for  large  inertia  (and  hence  reduced  motor 
damping) ,  the  larger  initial  overshoot  is  responsible  for  the  performance 
measures  being  outside  the  prescribed  bounds.  This  is  not  a  critical 
deficiency  since  the  initial  "attack"  of  the  response  represents  the 
cue  and  it  is  precisely  here,  rather  than  during  the  subsequent  washout, 
that  overshoot  can  be  tolerated.  On  the  other  hand  the  increase  in  IL 
for  the  low  amplitude  input  (0.25  w^)  and  inertia  of  2x10“  kg  m-2  is 
due  to  the  friction  effects  at  the  r  rirst  velocity  null  point.  Moreover, 
as  regards  friction,  two  observations  can  be  made:  firstly  for  large 
amplitude  reference  inputs  (such  as  w^)  the  velocity  lag  error  arising 
from  the  time  delay  between  input  and  response  at  w=0  is  sufficient  to 
overcome  static  friction.  Clearly  from  Figure  A1 (c)  this  is  not  the 
case  for  low  amplitude  reference  inputs  where  the  motor  velocity  remains 
zero  until  the  magnitude  of  the  velocity  error  is  sufficient  to  overcome 
stiction.  Secondly,  the  presence  of  friction  has  an  ameliorating  effect 
in  that  it  quenches  limit  cycles  that  would  normally  arise  from  encoder 
rounding  errors.  These  oscillations  are  evident  when  the  system  is 
simulated  without  friction  (Figure  A2(a))  and  are  reduced  in  amplitude 
when  the  encoder  wordlength  is  increased  from  10  to  12  bits  (Figure  A2(b)). 


5.  NONLINEAR  COMPENSATION 

The  results  of  simulation  indicate  that  the  effect  of  friction 
are  significant  for  low  amplitude  reference  inputs,  where  the  velocity 
error  at  zero  motor  velocity  is  insufficient  to  counter  static  friction. 

To  overcome  this, two  approaches  are  possible:  firstly  one  can  inject 
high  frequency  dither  at  the  motor  input  which  has  the  effect  of  "smoothing" 
the  nonlinearity  (10]  or  secondly,  apply  additional  nonlinear  compensation 
such  that  the  effective  loop  gain  of  the  system  increases  in  the  region 
of  zero  motor  velocity  for  low  amplitude  velocity  error.  In  the  present 
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case  direct  generation  of  dither  signal  by  the  digital  control  computer 
is  not  possible  due  to  the  relatively  low  sampling  rate,  and  hence  the 
alternative  amplitude  dependent  gain  approach  is  adopted  for  investigation 
The  development  below  should  be  taken  as  exploratory  rather  than 
definitive,  the  aim  being  to  demonstrate  the  flexibility  afforded  by 
computer  control  for  nonlinear  compensation. 

The  difficulty  of  detecting  the  null  velocity  point  of  the 
response, given  the  sampling  period  and  the  confutation  delay,  dictates 
that  the  nonlinear  control  effort  becomes  significant  within  some  region 
of  zero  motor  velocity,  for  low  velocity  errors  and  non-zero  reference 
inputs.  One  smooth  nonlinear  function  satisfying  these  criteria  is 


U  (k£t)  =  C  e  (kAt)  exp  (-a  e  (kAt)2  -  6  W(kAt)2), 
NL  W  r  w 

e  Aw  -  w 
w  =  r 


where  the  parameter  a  specifies  the  level  of  velocity  error  beyond  which 
linear  control  is  sufficient  to  overcome  static  friction,  g  specifies 
the  region  about  zero  motor  velocity  where  U  is  to  contribute  to  the 
control  effort  and  C  is  the  gain.  The  output,  u,  of  the  digital 
controller  can  now  be  taken  as 


u  (kAt)  =  uL  (kAt)  +  u^  (kAt) 


where  u^  is  the  output  of  the  linear  controller. 

The  choice  of  parameters  a,  g  and  C  are  determined  from  the 
simulated  response  of  the  closed  loop  system  to  the  velocity  input 
0.25  w  For  this  response  the  velocity  error  at  zero  motor  velocity 
is 


e 

w 


w=0 


1.5  rad/sec. 


Assuming  that  at  w-0  the  additional  torque  produced  by  the  nonlinear 
control  term  is  approximately  equal  to  the  static  friction  results  in 

k 

C(1.5)  exp(-Y  1.5*)  —  »  0.28  Nm  (s  stiction) . 

a 


Choosing 


exp(Y (1.5) a)  ■  0.8 


gives 


Y  -  0.1  and  C  -  1.56. 
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Further  the  region  about  w=0  where  the  nonlinear  control  is  to  be 
active  can  be  taken  as  (w|  <0.5  rad/sec  with 

1  ^  exp(-$wJ)  >  0.85  ;  [w|  <  0.5 


and  this  leads  to  - 

8  =  0.65. 


Results  of  simulation  for  these  parameter  values  are  shown 

in  Figure  A3,  where  it  can  be  noted  that  the  duration  of  zero  motor 

velocity  in  the  interval  about  t=1.8  sec  has  been  considerably  reduced 

compared  to  response  of  Figure  Al(c).  This  improvement  in  performance 

is  reflected  in  values  of  R  and  R.  as  recorded  in  Table  A2. 

w  w 

Little  analytical  machinery  exists  for  the  stability  of 
discrete-continuous  systems  with  multiple  nonlinearities.  The  discrete 
extensions  of  Popov's  Criterion  as  developed  in  til]  and  [12]  are  not 
applicable  to  the  present  case  as  they  are  limited  to  one  nonlinearity 
in  the  loop,  whilst  discrete  Liaponov  methods  ([3])  cannot  accommodate 
Coulomb  friction.  Nevertheless  it  can  be  agreed  that  the  closed  loop 
system  will  be  bounded-input  bounded-output  stable  since  for  large 
velocity  errors  and/or  large  motor  velocity  the  nonlinear  control  is 
insignificant  and  bounded  response  will  be  assured  by  the  linear 
analysis  of  Section  3.3. 


6.  CONCLUDING  REMARKS 


Frequency  domain  techniques  together  with  approximate  describing 
function  analysis  have  been  shown  to  offer  a  mathematically  tractable 
method  for  the  analysis  and  synthesis  of  digital  control  laws  for  a 
velocity  tracking  actuator  in  the  presence  of  load  variations  and  system 
nonlinearities.  A  significant  feature  of  the  control  law  is  the  inclusion 
of  a  nonlinear  gain  term  to  attenuate  the  deleterious  effects  of  stiction 
for  low  amplitude  inputs  without  adversely  affecting  the  satisfactory 
large  amplitude  response. 

The  principal  question  outstanding  relates  to  the  stability 
analysis  of  closed  loop  nonlinear  discrete  time  systems.  The  method 
presented  here,  based  on  continuous  time  describing  functions,  becomes 
cumbersome  and  of  limited  applicability  to  systems  with  multiple  non- 
linearities.  This  is  of  particular  importance  in  electromechanical 
systems  where  in  addition  to  common  nonlinear  effects  such  as  backlash 
and  saturation,  consideration  may  have  to  be  given  to  measurement  and 
control  quantization  effects  as  well  as  finite  computer  precision. 
Computationally  feasible  stability  techniques  based  on  discrete  (z-trans- 
form)  system  models  with  multiple  static  nonlinearities  need  to  be 
developed.  One  approach  warranting  investigation  is  the  method  of  Gossel 
(13]  based  on  generalized  superposition  principles  for  nonlinear  systems. 
This  could  offer  both  input/output  formulation  for  system  analysis  and 
a  mechanism  for  synthesizing  polynomial  type  nonlinear  control  laws  to 
achieve  magnitude  dependent  system  performance. 
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Performance  Measures  From  Simulation 
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Digital  control  strategies  for  a  nonlinear  motion  base  actuator 
are  considered  and  a  compound  linear/nonlinear  algorithm  is  derived 
for  velocity  tracking  under  a  range  of  load  conditions  for 
representative  motion  demands.  Discrete  frequency  domain  methods 
are  employed  to  synthesize  the  linear  components  of  the  control  law 
and  to  establish  the  stability  of  the  closed  loop  system  via  describ¬ 
ing  function  analysis.  The  nonlinear  compensation,  providing 
attenuation  of  Coulomb  friction  effects,  is  realised  through  a 
velocity  dependent  gain  term,  which  does  not  significantly  affect 
the  satisfactory  large  amplitude  system  response.  Overall  closed 
loop  system  performance  is  validated  through  computer  simulation. 
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